function [hu_fit] = hu_fit_wrapper_toolbox(patient);
    
    
    %|------------------------------------------------------------------------|
    %|Calculation of linear fit of HU with bellows voltage.
    %|  First slice registered images by slice
    %|  Then calculate linear fit each for slice to the deformed bellow
    %|  voltages;
    %------------------------------------------------------------------------
    %|      Dependancies;                                                     |
    %|                               |
    %------------------------------------------------------------------------
    %   This file is part of the
    %   5D-Novel4DCT Toolbox  ("Novel4DCT-Toolbox")
    %   DH Thomas, Ph.D
    %   University of California, Los Angeles
    %   Contact: mailto:dhthomas@mednet.ucla.edu
    %------------------------------------------------------------------------
    % $Author: DHThomas $	$Date: 2014/04/01 10:23:59 $	$Revision: 0.1 $
    
    if exist([patient.folder '/hu_fit.mat'])>0
        load([patient.folder '/hu_fit.mat'])
    else
        
        % Check if registered image *slices* exist, if not create them;
        patient.image_slice_folder = [patient.folder '/Deformed_Images_slices'];
        mkdir(patient.image_slice_folder);
        image_slice_dir = dir(patient.image_slice_folder);
        if size(image_slice_dir,1)-2<patient.dim(1);  %% Check if Registered Image Slices have already been created;
            create_deformed_img_slices_toolbox(patient)
        end
        
        hu_fit=hu_calibration_fitting_toolbox(patient);
%         hu_fit=fn_hu_calibration(folder, folder_static, ref, scans, static_mask, bellows_volt_drifted,directions,scan_times, run_par);
    end
